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ABSTRACT 

The separation of cells or particles from solution has traditionally been accomplished with centrifuges 
or by sedimentation; however, many particles have specific densities close to unity, making buoyancy-driven 
motion slow or negligible, but most cells and particles carry surface charges, making them ideal for elec- 
trophoretic separation. Both buoyancy-driven and electrophoretic separation may be influenced by hydro- 
dynamic interactions and aggregation of neighboring particles. Aggregation by electrophoresis was analyzed 
for two non-Brownian particles with different zeta potentials and thin double layers migrating through a 
viscous fluid. The results indicate that the initial rate of electrophoretically-driven aggregation may exceed 
that of buoyancy-driven aggregation, even under conditions in which buoyancy-driven relative motion of 
noninteracting particles is dominant. 

INTRODUCTION 

Electrophoresis is a common analytical technique for separating charged biological particles and macro- 
molecules. In an electric field, suspended particles will migrate at a velocity proportional to their surface 
charge and the applied field strength. Smoluchowski [1] showed that, for thin double layers, the elec- 
trophoretic velocity of an isolated particle is: 



where e is the dielectric constant, C is the zeta potential of the charged particle surface, Eoo is the applied 
electric field, p is the viscosity of the fluid, and the group in parenthesis is the electrophoretic mobility, 
p e , a physical constant of a given particle. Morrison [2] later showed that the electrophoretic mobility was 
independent of the size and shape of the particle. 

In a heterogeneous suspension, particles with different surface charges will have different electrophoretic 
velocities, leading to the possibility that they will collide and aggregate [3]. Such aggregation inhibits the 
ability for different subpopulations of cells or other biological particles to be separated and is the subject of 
this article. Particle aggregation has been investigated theoretically for buoyancy-driven motion by Davis [4] 
and Melik and Fogler [5] who used a trajectory analysis to determine the aggregation rates. A similar 
analysis for electrophoretic particle aggregation is performed in this work and compared with the results for 
gravity-induced aggregation. 

THEORETICAL DEVELOPMENT 

We restrict our attention to suspensions which are dilute, so that only two particles interact at one 
time. The particles are assumed to have negligible Brownian motion and the surrounding fluid is assumed 
to have negligible inertia. The zeta potentials are of moderate strength or less, and the electrical double 
layers are thin ( ko , » 1). These conditions are typically met for biological cells and other charged particles 
of 1-50 pm radius. 

To study coagulation and stability of dilute suspensions of spherical particles, a quantitative description 
of the relative motion between two spheres is needed. When the particles are far apart, this relative motion 
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is simply the difference between the separate velocities, given by: V° 2 = Vf — V^. For electrophoresis, 
V ® 2 = (1 — /?) Vf’°, where V B >° is defined by eq. (1), (3 is the ratio of zeta potentials, C 2 /C 1 - F° r gravity 
sedimentation, V ° 2 = (l — X 2 y) Vf’°, where V 6, ° is the Stokes sedimentation velocity [4]: 
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( 2 ) 


A p is the density difference between the particle and the fluid, a is the particle radius, g is the gravitational 
acceleration, A is the particle size ratio, a 2 /a i, and 7 is the reduced density ratio, Ap 2 /A/>i. 


As two particles approach each other, hydrodynamic interactions significantly affect the relative velocity, 
and this can be expressed through the use of mobility functions. Using the coordinate system shown in Fig. 1, 
we write [ 6 ]: 


Ux 2 ,r = - L(r ) cos e - ^ (, f A - N R f R ) 
Uiw = M(r) sin 0 , 


(3) 

(4) 


where Ui 2 ,r and Ui 2 t e are the components of the relative velocity directed parallel and perpendicular to the 
line of centers, respectively; the velocities are made dimensionless by Vj° 2 , and r is the distance between the 
centers of the two particles made dimensionless by the average particle radius. The mobility functions L(r) 
and M(r) describe relative motion due to the driving force. The mobility function G(r) describes relative 
motion due to interparticle forces, which are summarized later in this section. When the dimensionless 
interparticle separation distance, 6 = r — 2 , is large compared to unity, L(r), M(r), alnd G(r), approach 
unity, but, at near contact, L(r) and G(r ) approach zero, and M(r) tends to a finite limit. 

Figure 2 shows L(r) for electrophoretic and buoyancy motion. For electrophoresis, L(r) was obtained 
from the method of reflections [7] at large separations, a boundary collocation solution [ 8 ] at moderate 
separations, and a lubrication solution [9] at small separations. For buoyancy, L(r) is taken from the 
solutions summarized by Jeffery and Onishi [10]. The mobility functions for electrophoresis approach unity 
much more rapidly (with a leading-order correction proportional to 1 /r 3 ) than do those for buoyancy (which 
have a leading-order correction proportional to 1 /r) as the interparticle separation is increased; this is because 
electrophoresis is a force-free motion. 


When the interparticle separation distance is small, short-range interparticle forces will affect the relative 
particle motion. Electrostatic repulsion due to overlapping double layers can keep particles from aggregating, 
if strong enough. However, in most cases, van der Waals attractive forces will cause aggregation to occur. It 
will be assumed that these colloidal forces act only along the line of centers. The driving force for aggregation 
made dimensionless by the van der Waals force is described by the dimensionless parameter 


N a 




(5) 


where A is the Hamaker constant and a is the average particle radius. A dimensionless function describing 
the attractive van der Waals forces, taken from Gregory [11], is 


fA G^ 2 ) (5.325+1/) 


( 6 ) 


where 6 = h/a is the dimensionless particle separation distance, and v = A^/a is the dimensionless retar- 
dation parameter, with A^ « 0.1 fim being the London retardation wavelength. This expression is valid 
for 6 < 1, and, since van der Waals attractions are negligible at larger separations, eq. ( 6 ) will be used 
throughout. The dimensionless parameter describing the strength of the electrostatic repulsion relative to 
the van der Waals attraction is 

■ 


(7) 


380 



The dimensionless function describing the electrostatic repulsion force due to overlapping double layers, 
taken from Hogg et al. [12], is 


fn = 


f /3e- KaS - 

J . g— 2 na6 


( 8 ) 


This expression for the double layer repulsion is valid for na > 10 , Cl < 60 mV', and constant surface 
potential. 


RESULTS 


Colloidal suspensions tend to form aggregates of particles due to van der Waals attraction, but electro- 
static repulsion between two particles can lead to a stable suspension. To characterize the stability of the 
suspension, we examine the relative motion along the line of centers using eq, (3). A stable suspension is one 
in which the relative velocity along the line of centers, Ui 2 , r , becomes zero at finite separations for all values 
of —1 < cos 9 < 1. Moreover, the location at which Ui 2 , r = 0 is stable only if dUu^/dr < 0, indicating that 
a pair of particles will move apart if they become slightly closer than the stable separation distance. 

Melik and Fogler [5] have previously studied the stability of particles undergoing gravity sedimentation. 
We show here how their results can be extended to electrophoretic motion. In particular, if eq. (3) is written 
for gravity motion and electrophoretic motion, and set to zero relative velocity for each case, we have 

= *2 (9) 

along the neutral stability curves. Since L G JL E < 1 at small separations (see Fig. 2), this implies that 
electrophoretic motion becomes unstable at larger dimensionless driving forces than does buoyancy motion 
(Fig. 3). Moreover, since interparticle forces act at distances where near-contact hydrodynamics apply, 
L e /L g is independent of the small separation distance because L(r) is proportional to the gap between the 
particles for both types of motion [9]. 


Loewenberg and Davis [9] have analyzed electrophoretic near contact motion and compared the re- 
sults with buoyancy-driven near-contact motion. They find for small size ratios, A < 0.2, L E /L° can be 
approximated by 


L e _ 1.89 

L G “ A (1 — O.2O67) 


( 10 ) 


Then at small size ratios, L E /L G is proportional to A" 1 , which implies that electrophoresis is able to force 
particles together much easier than gravity is able to. In fact, L E fL G > 8 for all size ratios at 7 = 1 ; 
thus, electrophoresis is much less stable than buoyancy. Loewenberg and Davis [9] have postulated that 
the greater tendency for particles to aggregate when undergoing electrophoresis is due to electroosmotic 
convection which draws fluid out of the gap between two approaching particles. In contrast, the fluid must 
be squeezed out of the gap between two approaching particles in buoyancy-driven motion, which causes 
greater resistance to close approach. 


AGGREGATION RATE RESULTS 
Trajectory Analysis 

The pairwise aggregation rate is equal to the flux of particle pairs through an upstream collision cross- 
section [4]: 

J 12 = ni^ViWc , (11) 

where n< is the number density of particles of type i and <r c is the critical impact parameter at upstream 
infinity defined in Fig. 4 such that particle capture occurs only for <r < <r c . The collision efficiency is defined 
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as the actual collision rate given by eq. (11) divided by that in the absence of particle-particle interactions 
(<r c = ai + a 2 ): 



The critical impact parameter and the collision efficiency were determined from a trajectory analysis [4,5]. 
Dividing eq. (3) by eq. (4) yields 

dr__ (Hr) cos 6 + %g(f A -N R f R )) 

dO r M(r) sin 0 ’ ^ 

which was integrated backwards from a critical point at 6 = tt and the separation for which the driving force 
pulling the two particles apart is balanced by the interparticle colloidal forces pulling them together [4], 

Figures 5 and 6 show the collision efficiency as a function of the dimensionless driving force, N A , at 
different double-layer thicknesses and strengths of the electrostatic repulsion. Apparently, electrophoretic 
aggregation is approximately an order of magnitude more efficient than buoyancy-driven aggregation. The 
stability region (£12 = 0) for electrophoretic aggregation is shifted according to eq. (9), which is another 
manifestation that a given colloidal suspension is less stable to an applied electric field than to a gravitational 
driving force. As the driving force parameter increases, the collision efficiency at first increases as the 
electrophoretic or gravitational mechanisms pushing the particles together overcome the repulsive forces 
which hold them apart; a maximum is then reached, and further increases of the parameter N A lead to a 
decrease in £12 because the attractive van der Waals forces required for aggregation became relatively weak. 
By comparing Figs. 5 and 6, it is seen that adding electrolyte to decrease the double-layer thickness (increase 
ko) enhances the collision rate. This is because attractive van der Waals forces become stronger relative to 
repulsive electrostatic forces as the double-layer becomes thinner and shields the repulsive surface charge. 

CONCLUSIONS 

A suspension undergoing electrophoresis is less stable to aggregation than a comparable suspension 
undergoing gravity sedimentation. In addition, the collision efficiencies are much higher for electrophoretic 
aggregation than for gravity-induced aggregation. These results stem from the fact that electrophoretic 
motion is force-free and hydrodynamic interactions have a smaller effect on the relative motion than for 
gravity motion. 
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Fig. 1-Schematic of the relative 
motion of two charged spheres in 
an electric field. 


Fig. 2 -Relative mobility function along the line of 
centers for gravity motion (dashed line) and elec- 
trophoretic motion (solid line) with A = 2, 7 = 1, 
and 0 = 0.5. 


383 




Fig. 3-Stability plot showing sta- 
ble region, primary capture, and sec- 
ondary capture for gravity-induced 
aggregation (dashed lines) and elec- 
trophoretic aggregation (solid lines) 
with A = 0.5, 7 = 1, 0 1, v = 0.1, 

and Nr = 400. 



Fig. 5-The collision efficiency versus the ra- 
tio of driving force to attractive force for elec- 
trophoresis (thick lines) and gravity sedimen- 
tation (thin lines) for A = 2, 7 = 1,/? = 0.5, 
v = 0.044, Ka = 10, and Nr = 10 (solid lines), 
100 (dotted lines), and 1000 (dashed lines). 



Fig. 4-Schematic of the relative trajectories of two 
spheres for different impact parameters; <r c is the crit- 
ical impact parameter defined such that aggregation 
occurs for cr < <t c and not otherwise. 



Fig. 6-The collision efficiency versus the ra r 
tio of driving force to attractive force for elec- 
trophoresis (thick lines) and gravity sedimen- 
tation (thin lines) for A = 2, 7 = 1, /? = 0.5, 
v = 0.044, Ka = 100, and Nr = 10 (solid 
lines), 100 (dotted lines), and 1000 (dashed 
lines). 
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